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ABSTRACT 



Context. The Lyman alpha emission line of galaxies is an important tool for finding galaxies at high redshift, and thus probe the 
structure of the early universe. However, the resonance nature of the line and its sensitivity to dust and neutral gas is still not fully 
understood. 

Aims. We present measurements of the velocity, covering fraction and optical depth of neutral gas in front of two well known local 
blue compact galaxies that show Lyman alpha in emission: ESO 338-IG 04 and Haro 11. We thus test observationally the hypothesis 
that Lyman alpha can escape through neutral gas by being Doppler shifted out of resonance. 

Methods. We present integral field spectroscopy from the GIRAFFE/ Argus spectrograph at VLT/FLAMES in Paranal, Chile. The ex- 
cellent wavelength resolution allows us to accurately measure the velocity of the ionized and neutral gas through the Ha emission and 
Na D absoiption, which traces the ionized medium and cold interstellar gas, respectively. We also present independent measurements 
with the VLT/X-shooter spectrograph which confirm our results. 

Results. For ESO 338-IG04, we measure no significant shift of neutral gas. The best fit velocity is - 15± 16 km/s. For Haro 11, we see 
an outflow from knot B at 44±13 km/s and infalling gas towards knot C with 32±12 km/s. Based on the relative strength of the Na D 
absorption lines, we estimate low covering fractions of neutral gas (down to 10%) in all three cases. 

Conclusions. The Na D absorption likely occurs in dense clumps with higher column densities than where the bulk of the Ly a 
scattering takes place. Still, we find no strong correlation between outflowing neutral gas and a high Lyman alpha escape fraction. 
The Lyman alpha photons from these two galaxies are therefore likely escaping due to a low column density and/or covering fraction. 



Key words, galaxies: kinematics and dynamics - 
11 



galaxies: ISM - galaxies: starburst - galaxies: individual: ESO338-IG04 and Haro 



o : 

C<| 1. Introduction 



O 



x 



The Lyman alpha (Ly a) emission line was suggested as 
a probe for discovering high redshift galaxies already by 
iPartridge & Peeblesl (Il967l) . When absorbed by neutral hydro- 
gen gas, approximately two thirds of the ionizing photons from 
hot, massive stars are reprocessed into Ly a photons follow- 
ing case B recombination. The fraction of the bolometric flux 
contained in the Ly a line should be as high as 6-7 per cent 
in a young, star-forming stellar population. However, early sur- 
veys designed for findin g Lyman alpha emitting galaxies (LAEs) 
came up blank, (see e.g. lPritcheall994 and references therein). 

It was only with deeper and larger surveys (e.g. lCowie & Hul 

Il998l) that targeting Ly a became the successful method 
that it is today for finding galaxies at redshifts z > 
2. It has since then been widely used for mapping out 
the large-scale structure of the high re d shift universe (e.g. 
Rhoadsetal. 2000; Kudritzki et al. 2000; Malhotra & Rhoads 



2002 



Ouchi et al. 2003, 2005; Gawiser et al. 2006; Aiiki et al 



2006; Gronwa lleti alJ2007tlPirzkal et al.l2007tF inkelstein et al 
2008; Nilsso n et all l2009t lYumaetalJ l2010b lOuchietai 



* Based on observations made with ESO Telescopes at the Paranal 
Observatory under programme IDs 083.B-0470 and 60.A-9433 



20101 iGuaita et alJ[201 it lAdams et al.ll201 it iBlanc etai1l201 lb 
Shibuya et al. 2012 ) and even constraining the epoch of cosmic 
reionization (e.g.lMalhotra & Rhoadsl2004tlDiikstra et~aDl2007L 
lOno et alJlToiOtlJensen et al 112012b . 

As shown by lHaves et alj (|2010), up to 90 percent of star- 
forming galaxies in high redshift surveys emit too little Ly a to 
be detected by standard criteria. The first candidate to be blamed 
for this apparent discrepancy between theory and observation 
was absorption by dust, which is prominent in the ultraviolet. 
Early surveys in the local Universe hinted at an anti-correlation 
between the metallicity (which generally correlates with dust 
content) and Ly a luminosity (e.g. Meier & Terlevich 19811) . but 
it soon became clear that dust alon e could not explain the de- 
viation from recombination theory dGiavalisco et al.lll 996b . In 
particular one galaxy, I Zw 18, showe d a very low meta llicity 
combined with strong Ly a absorption (Kun th et alj|1994t) . 

The attention then turned to the resonant scattering of Ly a 
in H i, which had been explore d theoretically for a time (e.g. 
IOsterbrock|[l962t lAdamslll972"l) . If the path length of the Ly a 
photons is greatly increased in multiple scatterings on the way 
through the galaxy, even small amounts of dust can cause a large 
absorption (lNeufeldl l990). However, if the neutral gas is shifted 
in velocity, the Ly a photons are shifted out of resonance and can 
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escape from the galaxy more easily. iKunth et alJ (Il998l) showed 
from a sample of eight local gas-rich dwarf galaxies that the 
emission of Ly a in each case exhibited a P Cygni profile ac- 
companied by blueshifted low ionization state (LIS) absorption 
lines, suggesting that an outflow of neutral gas would allow the 
Ly a emission to escape towards us. M odern simulations (e.g. 
IVerhamme et al]l2006LlGarel et al.ll2012l Duval et al. 2012, sub- 
mitted) thus take into account a combination of the relative ve- 
locities of ionized and neutral hydrogen gas, as well as dust, and 
their respective morphologies, that all govern the escape of Ly a 
emission from galaxies. 

Although neutral gas is best traced through LIS absorp- 
tion lines in the UV, such a study at low redshift requires a 
space telescope due to ultraviolet absorption in the atmosphere. 
The choice is then either HST/STIS which has a low sensitiv- 
ity or HST/COS with limited spatial information. However, the 
emission and absorption of Ly a may va ry over small scales 
(M as-He sse et a U mM. lOstlin et al] l2009) and a detailed study 
requires some degree of spatial resolution. It is therefore use- 
ful to identify an alternative set of longer wavelength absorp- 
tion features that still form in the neutral ISM, in order to facil- 
itate the use of contiguous integral field spectrographs attached 
to large-aperture telescopes on the ground. Neutral gas motions 
in star-forming galaxies h ave been studied in the past using, 
e.g. Mg n 2796,2803 (e.g. IChurchill & Vogj|200ll: i Mshar et al l 
120071: iMartin & Bouchel |2009| iNestor et al.1 l201ll) and Na i 
5889. 5895 (e.g.lHeckman et al.ll2000l:lRupke et alJl2002l:lMartid 



2005| IChen et al.ll2010ir although never before have such kine- 
matic tracers been used with reference to Ly a. 

In this paper, we present the first spatially resolved study of 
absorption from Na i in the cold interstellar gas in two nearby 
Ly a emitting blue compact galaxies. Specifically, we target 
the sodium resonance absorption doublet (AA 5889.95,5895.92), 
which we will refer to as Na D in this paper. The galaxies are 
Haro 1 1 (ESO 350-IG38) and ESO 338-IG04 (Tololo 1924-416) 
which are both objects of intense study due to their Ly a emis- 
sion and the possible similarities to LAEs and Lyman Break 
Galaxies (LBGs) at high redshift (see Section|2]i. 

Na D is generally stronger than any other resonance line in 
the optical (such as K i or Ca n, which are strongly depleted in 
diffuse, low-velocity clouds, see e.g. ISpitzeri ([1968)). With an 
ionizing potential of 5.14 eV it is a good tracer of neutral hydro- 
gen gas. By simultaneously studying the Ha emission line, we 
can accurately determine the location, extent and velocity dis- 
tribution of the ionized hydrogen in these galaxies. We are thus 
able to measure the relative velocities between the ionized and 
neutral hydrogen gas, which is precisely the Doppler shift rele- 
vant for Ly a transmission. With integral field unit spectra from 
Argus at VLT/FLAMES we achieve a high spectral resolution 
combined with spatial information. We include an independent 
measurement from VLT/X-Shooter to confirm these results, and 
to estimate the stellar contamination of the spectra. 

The paper proceeds as follows; in Section 2 we describe pre- 
vious studies of the two galaxies in our sample in more detail. 
In Section 3 we describe our data and the reduction steps. In 
Section 4 we present our results. In section 5 we discuss pos- 
sible complications from stellar contamination. In Section 6 we 
discuss our results and in Section 7 we leave our concluding re- 
marks. 

2. Sample 

The two galaxies in our study have already been studied exten- 
sively for their interesting properties, but the neutral gas veloc- 



ity as estimated from the sodium doublet has not been measured 
with high accuracy and spatial resoluti on until now . Anal yzing 
VLT/UVES spectra of ESO 338-IG04. lo"stlin et alJ (12001 esti- 
mated an outflow velocity ~ 20 km/s from Na D, but emphasized 
that with their low signal to noise they could not estimate the rel- 
ative ISM versus stellar contribution. 

An HST UV i maging progr am includin g the two ga laxies 
(IKunth etalJliOOlt lHaves et alJl2005l 12001 : lOstlin eUu1l2009l) 
gives us information on the Ly a and UV continuum morpholo- 
gies. 

Haro 1 1 is a blue compact galaxy at z = 0.02 with three large 
condensat ions or "knots". T he knots are traditionally called A, 
B and C dVader et al.ll 19931) as shown in the left panel of Figure 
[1] Knot C appears to be the brightest knot in the ultraviolet and 
it also shows Ly a in emission, while knot B is brightest in Ha, 
where Ly a is instead absorbed. These two knots are therefore 
of particular interest as we try to explain why Ly a is seen from 
one region but is absorbed in the other. Some evidence suggests 
that Haro 1 1 might b e the result of a merger of dwarf galax- 
ies dOstlin et al. 2001), particularly due to the irregular appear- 
ance, the high relative velocities of several hundred km/s, the 
broad emission lines, and the presence of a tidal arm structure 
with a high redshifted velocity relative to the rest of the system 
(Ostlin et al. 2013, in preparation). The H a line width within a 
star forming knot is as high as ~ 270 km/s (FWHM), and shows 
strong multi-component features. 

IKunth et alJ dl998h present HST/GHRS spectra of Haro 1 1 
around Ly a and the interstellar O i (A 1302 A) and Si n (/H304A) 
lines, estimating an outflow velocity of the neutral medium of 
circa 60 km/s. Both absorption lines are very broad, indicating 
multiple components along the line of sight spanning roughly 
200 km/s in velocity. The Ly a line profile has a strong underly- 
ing absorption, extending more than 1500 km/s on the blue side 
of the emission. The line does not show a clear P Cygni pro- 
file, and the underlying absorption seems to extend also to the 
red side. Evidence for several blueshifted Ly a components was 
found, indicating multiple absorbing gas clouds. Unfortunately, 
it is not clear exactly where in the gala xy the l'/7 x 1"7 GHRS 
aperture was pointed (see discussion in Haves et al. 20071). 

ESO 338-IG04 (or just ESO 338 as we might refer to it in 
this paper) has slightly narrower emission lines than Haro 1 1 
(FWHM ~ 180 km/s) and the Ly a emission is concentrated to a 
bright central region. ESO 338 is riddled wit h "super star clus- 
ters"; small knots of intense star formation (lOstlin et al.l f 1998. 
120031 120071) . The ongoing starburst is about 40 Myr old and 
was likely triggered by a merger with a small galaxy or from 
an interaction with remaining debris from a previous encounter 
with the companion galaxy ESO 338-IG04b (lOstlin et alj|2001t 
Cannon et al. 2004). When we discuss ESO 338 below, we 
mainly focus on the bright central region, which we refer to as 
knot A (lHaves et al. 2005[ cf. Figure [0 (also k nown as cluster # 
23 in the 'inner' sample of lOstlin et all (1 19981) ). 

lHaves et"ail d2005l) present a detailed description of a 
HST/STIS long slit spectrum across ESO 338. The spectrum 
shows diffuse Ly a emission in several regions along the slit, but 
over knot A the Ly a emission is weak and shows only a hint of 
blueshifted absorption. However, the low resolution of the spec- 
trum makes a detailed kinematical study difficult, and the emis- 
sion line lies very close to the strong geocoronal Ly a line. The 
slit is also only U.'2 wide, and was not optimized to be centered 
specifically over knot A. It is therefore not clear how well this 
spectrum samples the Ly a profile in this region. An outflow ve- 
locity from this STIS spectrum was estimated by Schwartz et al.l 
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Fig. 1. Slit positions, field of view (FOV) and approximate aperture positions for Haro 1 1 (left) and ESO 338-IG04 (right), overlayed 
on HST/ACS F550M continuum images (i.e. close in wavelength to the Na D feature). The X-Shooter slits are 0.9 (or 1.0 for the 
UVB-arm, not marked in the figure) arcsec wide and 1 1 arcsec long and marked with solid lines. The FOV for the LR6 wavelength 
range (which includes Ha) is marked with dash-dotted lines while the shorter LR5 range (including Na D) is marked with dotted 
lines. The Argus binned "apertures" are 1.56 arcsec wide squares marked with dashed lines. The VLT/FLAMES FOV shown here 
shows the area for which we have full exposure after taking into account the dithering pattern. 



(2006) from low ionization state absorption lines in the UV to be 
47 + 70 km/s. 

In this paper we look mainly at the strong knots in the respec- 
tive galaxies; see Figure Q] Even though our field of view with 
VLT/FLAMES is larger, the Na D signature is very weak and the 
areas marked are the only regions where we can successfully fit 
the Na D lines, even after spatially binning our spectra. 

3. Observations and Data reduction 

3.1. Integral Field Spectroscopy with VLT/FLAMES 

The integral field spectroscopy observations (ESO ID 083. B- 
0470(A)) were performed between June 9 and July 6 2009 with 
the Argus integral field unit of the GIRAFFE spectrograph at 
the FLAMES instrument at VLT/UT2. Argus consists of 14x22 
lenslets arranged in a rectangular grid, corresponding to a field 
of view of approximately 7x11 arcseconds. Each lenslet is then 
connected to a fiber which leads the light to a spectrograph. One 
spectrum for each lenslet is thus obtained. 

We observed the two galaxies Haro 1 1 and ESO 338-IG04, 
both with two "low resolution" (H = 11800 - 13700) gratings. 
The LR6 grating (which covers A 6438-7184 A) was used to 
observe the Ha emission line (A 6562.82 A), and the LR5 grating 
(A 5741-6524 A) at the Na D doublet (AA 5889.95,5895.92 A). 

The exposure time was 195s x 8 = 1560s for the LR6 grating 
and 895s x 12 + 915s x 2 = 12570s for LR5 for both galax- 
ies. The seeing was 075 - 078 when using the LR5 grating and 
for LR6 it was ~ 078 for Haro 1 1 and 079 - 173 for ESO 338. 
The basic bias, flat, wavelength calibrations and sky subtractions 
were made with recipes from the common pipeline library (CPL) 
package version 5.2.0 in esorex. 

The bias frames were median combined with the 
gimasterbias recipe, producing master bias frames for 
each observing night. The gimasterf lat recipe was then used 
for calibrating the fiber positioning onto the CCD. The flat 
frames used for each observing block were always acquired 
during the same night, using the Nasmyth screen (which gives 
better illumination in Argus mode). The giwavecalibration 



recipe subsequently performs the wavelength calibration by 
mapping known wavelengths from a calibration lamp. This 
recipe requires an existing dispersion solution as an initial 
guess, which is then refined. We checked the convergence of 
this procedure by giving the refined dispersion solution as a new 
initial guess and ensuring that the two refined solutions were 
identical for each calibration set. Finally, the giscience recipe 
applies all calibrations to the data and creates a data cube with 
two spatial and one spectral dimension. 

Each science spectrum was reduced individually. As a sanity 
check, we verify the sky subtraction by measuring the continuum 
flux in regions far from the bright condensations, which we find 
to be consistent with zero. 

The subsequent reduction steps were made with the PyFITS 
module version 2.3.1 for Pythory. Since the science data were 
dithered (the telescope made small movements of roughly 0.5 to 
1 arcsec between each science exposure), we created an empty 
x-y-/i-cube with spatial dimensions matching those of the region 
corresponding to the region on the sky for which we had con- 
tributions from each science exposure. The spectral data were 
then combined onto this grid using a weighted average intensity 
for each integral field unit pixel (or "spaxel"), calculated from 
the uncertainty values obtained from the pipeline. The dithering 
shifts were rounded off to whole spaxels before combining. This 
may thus introduce a small astrometric shift of up to 0725 in an 
individual exposure. Our final spatial resolution is therefore ex- 
pected to be similar to our worst seeing of about one arcsec. The 
final data product has the original 0752/lenslet spatial sampling 
resolution and 0.2 A spectral sampling resolution. However, we 
always consider binned spaxels in a 3x3 configuration in this pa- 
per, corresponding to a bin size of 1756. 

Line profiles where fitted using the optimize. leastsq routine 
from the Scip>Q version 0.7.0 Python package, which uses a 
Levenberg-Marquardt least-squares minimization technique. 



1 PyFITS is a product of the Space Telescope Science Institute, which 
is operated by AURA for NASA 

2 http://www.scipy.org 
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Table 1. X-Shooter observing specifications Table 2. Na D properties and dust extinction 



Parameter 


Haro 11B 


Haro 11C 


ESO 338-IG04 


Airmass 


1.016 


1.012 


1.430 


Seeing 


0T7 


a/7 - 1'/3 


a/ 6 


Exposure time 


800 s 


680 s 


600 s 



3.2. VLT/X-Shooter spectra 

The X-Shooter data were obtained as part of the first science 
verification for the instrument between August 10 and 11 2009 
(ESO ID 60.A-9433(A)). For the VIS arm (A ~ 5500 - 10200), a 
slit of a/9 x 1 1" was used. For the UVB arm (A ~ 3000 - 5600), 
the slit was 1" x 11". The resolving power % is 8800 and 5100 
for the VIS and UVB arms, respectively. 

The spectra were reduced with the X-Shooter pipeline v. 
1.3.7 using esorex v. 3.9.0. Standard settings were used in the 
physical model mode, using the xsh_scired_slit_nod recipe 
to perform reduction, sky subtraction and extraction of the sci- 
ence frames. 

These data are also analyzed bv lGuseva et al.l (1201 2l) . but we 
have included different reduction steps and we do not discuss the 
data from the NIR arm. With the large wavelength range of these 
data, we can examine the effect of the underlying photospheric 
Balmer absorption on the Ha emission feature. It is strongest in 
ESO 338, but negligible for our discussion in all cases. 

4. Results 

The Na D feature is very weak in these galaxies and we were 
therefore forced to bin our FLAMES spectra in the spatial di- 
mension. Even in this case, only the very strongest continuum 
features show Na D absorption. This reduces our discussion of 
the larger fields of view of FLAMES to the selected binned aper- 
tures shown in Figure [TJ Except for knot A in Haro 11, these 
are the only regions where we can see a clear Na D absorption 
profile. Nevertheless, these apertures still allow us to perform a 
sufficiently detailed analysis in some of the most interesting Ly 
a emitting and absorbing regions. 

The most important spectral features from the FLAMES 
binned apertures are shown in Figure |2] The spectra are not flux 
calibrated and are instead shown normalized to the continuum. 
The continuum was in all cases fitted with a low-order poly- 
nomial across the entire observed range, which yielded a good 
fit. The Ha and He i (5875.64 A) lines both show strong multi- 
component features, as is evident from the deviation from a sim- 
ple gaussian. Note that these features are very similar in both 
emission lines. We verify that the asymmetry is not an instru- 
mental feature by examining the shape of skylines and of cal- 
ibration lamp frames. The measured velocities of the different 
spectral components are listed in Table [3] 

We have also attempted multi-component fits to the Ho- emis- 
sion lines (not shown). In all cases, we identify two strong com- 
ponents in the Ha lines: One narrow "main" component and 
one broader but weaker component. In the case of Haro 1 1 B 
and ESO 338, the velocity of the broad component is consistent 
with the velocity of the main component, within the uncertain- 
ties. For Haro 1 1 A the narrow (broad) component has a velocity 
of 6242 ± 7 (6189 ± 10), and it is clear that the profile is very 
complex, but we have no good Na D velocity to compare it to. 
In the case of Haro 1 1 C the weak, broad component is centered 
roughly at the same velocity as the Na D profiles. In this case, 
the contribution to the Ly a flux from this component is there- 



Parameter 


Haro 11B 


Haro 11C 


ESO 338-IG04 


Na D doublet EW (A) 


-0.36 ± 0.06 


-0.42 ± 0.06 


-0.15 ±0.04 


Na D 5890/5896 ratio 


1.15 ±0.22 


1.14 ± 0.18 


0.94 ± 0.47 


Na D optical depth 


> 2 


> 2 


> 1.6 


Na D mean shift (km/s) 


-44 ± 13 


32 ± 12 


-15 ± 16 


Nebular E(B-V) 


0.42" 


0.48" 


< 0.1* 



Notes. References: al lHaves et aP d2007l) . b) Berevall lTOstlirJ 12002) 



fore not only expected to be weaker intrinsically but should also 
suffer a larger resonance effect. In each case, these secondary 
components do not change the velocity of the main component 
by more than a few kilometers per second and for simplicity we 
have chosen not to include them in the further discussion of this 
paper. Our goal is not to perfectly reproduce the line shapes but 
to analyze the velocities of the main constituents of the galaxies. 
We note that the componen ts we measure in H aro 1 1 agree well 
with the recent results from Jame s et al. l (l20H . 

We also attempted a simple wavelength centroid fit, numeri- 
cally defined as "LAfxI^fx across the continuum subtracted emis- 
sion lines. The result is perfectly consistent with the numbers we 
use in this paper, with a maximum shift in velocity of 5 km/s. 

The X-Shooter spectra are shown for comparison in Figure |3] 
(see also Section|5]l. Note that single gaussian fits appear to re- 
produce the line shapes better for these data, but this is mainly 
due to a lower spectral resolution and signal to noise. 

We see the velocity of emission lines from ionized gas are the 
same within each region we measure, independent of the species 
used. In the same region, the velocities as measured by the emis- 
sion lines of Ha (A 6562.8 A), He i (A 5875.6 A), [O i] (AA 
6300.2, 6363.9 A) [N n] {AA 6548.1, 6583.6 A) and [S n] (AA 
6716.5, 6730.7 A) are all the same, within uncertainties. 

In general, the Na D features are weak in both galaxies. For 
Haro 11, we measure an equivalent width for the doublet of 
EW^d = -0.42 ± 0.06A for knot C, -0.36 + 0.06A for knot B, 
and -0.15 + 0.04A for ESO 338. For Haro 1 1 A we estimate an 
upper limit of -0.15A. These values are estimated purely from 
the FLAMES spectra. The ratio between the EW of the separate 
lines of the doublet appears relatively close to unity (cf. Tabled, 
which indicates a high optical depth in the individual clouds that 
fall on sight lines where absorption is occurring. Note, however, 
that this is not the same as the average optical depth along all 
sight lines in the aperture, which is clearly smaller. The Na D line 
ratio in the optically thin case is 2:1. We present in Table |2]the 
estimated optical depths at the line center of Na D, deriv ed from 
the li ne ratio (using the conversion in Table 2.1 from (jSpitzer 
1 19681) ). which are all found to be optically thick. In this regime, 
the covering fraction can be easily estimated as Cf = 1 - /5890 
where /5890 is the residual intensity in the blue Na D line. Note 
however, that the residual intensity may be low in narrow, unre- 
solved components in our spectrum and that we therefore mea- 
sure an artificially stronger intensity. We will come back to this 
point later in the discussion. With this simple assumption how- 
ever, we estimate as a lower limit the covering fraction of Na D 
to be ~ 10 % for Haro 1 1 knots B and C, and roughly 5 % for 
ESO 338. 

For Haro 11, we detect no measurable Na D absorption to- 
ward knot A. This is probably due to the low continuum bright- 
ness in this region which is necessary for observing absorption 
lines, combined with a low covering fraction of neutral gas. For 
knot B, we are able to measure the Na D velocity and find a 
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Fig. 2. Hot, He i (5875.6) and Na D line profiles from the VLT/FLAMES data. The blue gaussians are a simple fit to the emission 
or absorption lines. In the case of Haro 1 1 knot A, the best fit is shown in this figure but is never used. The vertical dashed lines 
show the velocity given by the Ha gaussian fit, and the vertical dotted lines show the Na D velocity. The y-axis is normalized to the 
continuum. Also shown inset in the bottom row are errorbars, representing the standard deviation in the continuum from 5900 to 
5950 A which shows no emission or absorption features. Since the spectra are normalized to the continuum, these errorbars reflect 
the strength of the continuum in each region. 



moderate blueshifted velocity (compared to the ionized gas) of 
44+13 km/s. In knot C the absorption is stronger and we find a 
redshifted velocity of 32±12 km/s. 

In ESO 338, we are only able to measure the Na D velocity 
in the brightest knot (knot A), were both the continuum and the 
Ly a emission is strongest. We find no evidence f or a shift in 
the ve locity, with the best fit giving -15+16 km/s. lOstlin et alJ 
J2007h estimated a shift of circa 20 km/s from UVES spectra 
with a lower signal to noise. Our estimate of th e stellar versus 
nebular velocity agrees well with the analysis of Cumming et alJ 
(2008), where the stellar component (as measured by the Ca n 
triplet (AA 8498, 8542, 8662 A) shows the same velocity as the 
nebular gas, with a value of 2860+4 km/s. 



5. Possible stellar contamination of Na D 

A challenge with using Na D to measure neutral gas flows is that 
the stars in the galaxy may also exhibit the absorption lines. In 
this section we demonstrate that the Na D we measure is primar- 
ily of interstellar origin. 

Na D is prominent in spectra of cool stars: K- and M-type 
giants and supergiants. There are several ways of disentangling 
the stellar and interstellar sodium. One of the more common and 
robust methods make use of the Mg i b-band triplet (AA 5167.32, 
5172.68, 5183.60 A). It has a very similar ionization potential 
as Na, and the two elements are formed under similar conditions 
in stars. However, the triplet forms between excited states and is 



not found in the cold, neutral medium. Therefore, any detection 
of Mg i b must come from stellar atmospheres. 

Thanks to their similar origins, the equivalent widths of 
the Mg b triplet and the stellar Na D doublet appear to 
be strongly correlated. However, any ground-based spectra of 
nearby stars are always contaminated by telluric Na D in 
emission, which m akes the relative ratio difficult to measure. 
iRupke et alJ J2002) compile measure ments of Galactic glob ular 
clusters and mostly nonactive galaxies dBica & Alloinll986h and 
nuclei of nearby galaxies dHeckman et al.l H980). They find the 
equivalent widths of the purely stellar features to be correlated as 
EWnhd ~ 0-5 EW Me b with a possible intrinsic s catter of > 25%. 
iHeckman et alJ (l2000h combine the d ata from Heckman et alJ 
(Il980h with stellar library data from iJacobv et al.l dl 984-1) and 
estimate EW^d ~ 0.75 EWm s b> which s eems to agree with 
their o wn K giant spectra to within 0. 10 dex.lSchwartz & M artin 
(2004j) also make a fit to the data bv IJacobv et all d 19841) to 
find EW NiD = 0.40 ± 0.05 EW Ug b. iMartinl d2005"h analyze 
Keck II/ESI spectra of A, F, G and K dwarfs and giants and 
conclude EW^d — 1/3 EWy^ % b and that the stellar contri- 
bution w as consistently le ss than 10 % in their study of 18 
ULIRGs. ISato et all d2009h investigate Na D absorption in 493 
spectra from the AEGIS survey, and find £WNaD = 0.40 EWMe, b 
well describes the p urely stellar bound ary in their sample (see 
their Fig. 1). Finally. IChen et al.l (l2010h used Sloan Digital Sky 
Survey (SDSS) spectra from young disc galaxies and estimate 
that an average ~ 80% of the Na D absorption arises in stellar 
atmospheres. Their estimate of the stellar contribution is based 
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ESO 338 



Table 3. Best fit observed velocities, in km/s. 
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Fig. 3. Ha, He i (5875.6), Na D, K i and Mg i line profiles from 
the VLT/X-Shooter data. The blue gaussians are a simple fit to 
the emission or absorption lines. The vertical lines show the sys- 
temic velocity as given by the Ha gaussian fit. The y-axis is nor- 
malized to the continuum. 
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Notes. Mg i and K i are measured only from VLT/X-Shooter and the 
uncertainty is therefore larger. The other lines were measured in both 
VLT/FLAMES and VLT/X-Shooter and this represents the weighted 
average. 



lower resolution of the X-Shooter data, the equivalent widths are 
not as easily determined as for the Argus data. However, we can 
still use our spectra in Figure [3]to see that it is only in Haro 1 1 
C that we can clearly distinguish the Mg b feature, although it is 
very weak. We thus let the upper limit from Haro 1 1 C represent 
the upper limit on the stellar contamination for all our results. If 
the EW of Mg b would be larger than this limit in the other re- 
gions this would be readily seen in Figure[3] We estimate the EW 
of this feature to be roughly -0.2 A, i.e. about half the equivalent 
width of Na D. 

Based on the discussion above we thus estimate the stellar 
contribution to be on the order of 25-30%, but we note that it may 
be as high as 50% in the extreme case. A significant contribution 
to the Na D profile at the stellar velocity (which is consistent 
with the H-alpha velocity in our galaxies) would decrease the 
average offset velocities that we measure. The absolute value of 
our Na D velocities may therefore be somewhat lower than the 
purely nebular Na D velocity. 



on fitting and subtracting a stellar population synthesis model to 
the continuum and absorption features, and they emphasize that 
their spectral model sometimes shows stronger Na D absorption 
than the actual data. They conclude that the most likely expla- 
nation is that Na D is sometimes seen in emission in their data, 
predominantly from face-on galaxies with low dust attenuation, 
and that a stronger stellar contamination is expected in the spec- 
tra of normal star-forming galaxies as opposed to that in young 
stars comparable to the dominant stellar population in our two 
galaxies here. 

We note that Na D is a resonance absorption feature, and 
therefore absorbed photons will be re-emitted in a random di- 
rection. This can create a diffuse , low surface brightness emis- 
sion component. iProchaska et all (1201 ll) model the effect of re- 
emitted light and conclude that in the extreme case that all of 
this light eventually escapes and is caught within the aperture, it 
can significantly reduce the observed equivalent width. We ex- 
pect this effect to be very small in our analysis, however, since 
we are looking only at narrow slits and small effective apertures 
(1.5 arcsec corresponds to roughly 0.6 kpc in Haro 1 1 and even 
less in ESO 338). 

An alternative to measuring outflows with Na D is the K i 
doublet (AA 7664.91, 7698.97 A). With similar properties and 
a similar ionization potential of 4.34 eV (compared to 5.14 eV 
for Na D), th ese lines are expec ted to behave much in the same 
way (see e.g. Kemp et al. 2002). They are however very seldom 
used, because they are weaker than Na D due to a lower relative 
abundance but also because they often overlap with atmospheric 
O2 absorption bands. 

Here, we use our X-shooter data to estimate the stellar con- 
tamination. Due to the comparatively low signal-to-noise and 



6. Discussion 

We have measured outflow/inflow velocities and covering frac- 
tions of the neutral ISM in front of three bright star-forming 
condensations, two of which show Ly a emission (ESO 338 A, 
and Haro 1 1 C) and one shows Ly a in absorption (Haro 1 1 B). 
Further to the information measured in this programme, we also 
compile measurements of the nebular d ust attenuation in th ese 
regions fr om our previous investig ations: lHaves etail (120071) for 
Haro 11, Berg vall & Ostlinl d2002l) for ESO 338, see Tableg] 



6.1. Comments on the individual regions 

Emission - ESO 338 knot A For ESO338-IG04, we find that 
the velocity of the neutral gas (as estimated from the Na D ab- 
sorption) is only very slightly blueshifted (or even static) com- 
pared to the ionized gas. This would rule out the outflow sce- 
nario as an explanation for the observed direct Ly a escape 
from this region. However, we note that Ha from the central 
regio n (knot A) of ESQ 338 is dominated by a large H 11 shell 
(e.g. Bergvall & Ostlinl 120021) . The UVES spectra analyzed by 
Ostl in et alJ d2007l) show a multi-component feature in the [O 
in] (/15007A) and Ha lines towards knot A, which they interpret 
as an expanding bubble at ~ 40 km/s. Presumably this shell does 
consists of outflowing ionized gas from stellar winds and super- 
novae, yet we measure a low outflow velocity in H 1. If there is 
indeed an expanding ionized bubble in ESO 338, it is possible 
that the Ha feature is dominated by emission on the side of the 
bubble facing our way. This would cause us to measure a more 
blueshifted velocity for the ionized gas and reduce our inferred 
neutral gas velocity. However, the bubble seems to be optically 
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thin, which would mean that the effect is very small, and also the 
stellar component i n knot A appears to have the same velocity 
as the ionized gas dCumming et al.ll2008l) . 

Emission - Haro 1 1 knot C For Haro 1 1, we find a small red- 
shifted velocity of Na D in front of knot C, indicating infall of 
cold gas. HST imaging shows Ly a emission from this region, 
which could partially be explained by this velocity difference. 
However, it is in the opposite direction to an outflow, and the 
velocity difference is not particularly strong. 

There is evidence for a Ly a / Ha ratio high er than the rec om- 
bination case B value in both ESO 338 A dOsflin et al.l2009l) an d 
Haro 1 1 C (if extinction is taken into account: [Atek et all l20Q8). 
This could be due to Ly a photons actually suffering less at- 
tenuation while scattering on the surface of cool dusty clumps 
than the continuum photons (Neufeld 1991), alth ough this sce- 
nario seems unlikely based on recent simulations (lLaursen et alj 
1201 2l Duval et al. 2012, submitted). It is also possible that an 
extinction correction based on a clumpy dust distribution model 
rather than a uniform dust screen w ould give a result co nsistent 
with recombination values (see e.g.jScarlata et al. 20Q9|). 



neutral gas, consistent with the strong Ly a absorption on the 
blue side of the line. The resolution of GHRS {ft ~ 16000) is not 
much higher than that of our Argus Na D spectra (H = 1 1800), 
and would likely not explain this difference. We note also that 
other Na D absorption spectr a in the literature often show con - 
siderably stronger lines (e.g. lHeckman et al.l2000llMartinl2 005). 
Our conclusion is that the Na D covering fraction is indeed low 
in these regions. The difference from the previous measurement 
might in part be explained by the different apertures pointing 
at different sightlines in an inhomogenous ISM. The difference 
may also come from the lower photoionization threshold of Na 
D of 5. 14 eV. It is thus possible that the column density of hydro- 
gen is high enough overall to absorb Ly a but low enough that we 
detect Na D only in the densest regions. The interstellar medium 
is known to often be very patchy, and it would not be surprising 
in these two galaxies with their turbulent past s. Since Na D has 
to be shielded from ionizing radiation by dust (IChen et alJ l2010; 
iMurrav et al .1120071) . it likely only exists in the densest, coolest 
clumps. 

For the optical depths that we find from the Na D line ratio, 
we can put a lower limit on the column density of the gas where 
the absorption arises. We use the relation 



Absorption - Haro 1 1 knot B For knot B, we find a moderate 
blueshifted velocity. Combined with a low covering fraction ~ 
10 %, Ly a should escape more easily from knot B. However, 
our HST UV imaging shows that Ly a is strongly absorbed in 
this region. 

There is some debate as to whether the dust extinction is 
higher in knot B than in knot C. Judging from HST images, 
there do seem to be more dust clouds near and around knot B 
and based on the X-Shooter d ata, knot C seems to have a lower 
extinction dGuseva et al.l2012h . However, the Ha and H/3 images 
from which we derive the E(B-V) values presented in Table [2] 
give approximately the same level o f extinction in both knots 
(iHaves et alJ 120071 lAtek et al]l2008l) . If the medium is highly 
clumpy it is likely that the extinction is varying, even within 
these small regions. Different results are thus likely due to differ- 
ent aperture sizes and positions, as well as different techniques 
for measuring the extinction probing varying optical depths. 

Absorption - Haro 1 1 knot A The signal-to-noise and the 
equivalent width of Na D are too small for us to safely attempt 
a measurement toward knot A. This is partly because the optical 
continuum emission near Na D is weaker than in the other knots, 
but it still indicates a low covering fraction of neutral gas. Knot 
A exhibits both H a and UV continuum emission from ionized 
regions, but no significant Ly a radiation appears to escape from 
them. 



6.2. Interpretation of the covering fractions, optical depths 
and velocities 

In the three regions that we have identified for study in this arti- 
cle, we have compiled measurements of H i covering fractions, 
kinematics, and dust contents. 

In both knot A and B in Haro 11, absorption of Ly a is seen 
despite a low covering fraction. However, the covering fractions 
that we estimate from the Na D profiles likely only serves as 
lower limits, since there may be very narrow, unresolved com- 
ponents with lower residual intensity in our s pectra. We note that 
the O i and Si n lines in the GHRS spectra of lKunth etal] d!998h 
are considerably stronger, indicating a high covering fraction of 



N(Na i) = 



Tib 



1.497 x 10- 15 Ai/i 



(1) 



from lSpitzeri (Il978l) . where t\ is the central optical depth, b 
is the Doppler parameter (in km/s), A\ is the vacuum wavelength 
in A and f\ is the oscillator strength = 0.3180, all for the weaker 
(red) line. With a Doppler parameter b = FWHM/(2 Vln2) of 

> 60 km/s for the Na D in ESO 338, a lower limit on the Na i 
column density is roughly N(Na i) > 4 x 10 13 crrT 2 . Converting 
this to a hydrogen column density of course depends heav- 
ily on the Na / H abundance ratio, but using ve ry conserva- 
tive es t imates based on the conversions given in Ru pke et al.l 
d2002l) : IMurrav et al.l (120071) this corresponds to at least N(H i) 

> 10 20 crrT 2 . Ly a becomes optically thick and resonantly scat- 
ters on neutral hydrogen already at H i column densities around 
10 13 ~ 14 crrT 2 . It would take an enormous departure from these 
assumptions to make the clouds optically thin. In these dense 
clouds, the Ly a photons would see upwards of a million op- 
tical depths, and the clouds would be self-shielding. Still, the 
covering fraction of these dense clouds where we see Na D is 
only ~ 0. 1 . Our interpretation is that the ISM in these galaxies 
is inhomogeneous and likely consists of column densities be- 
tween these extremes. Ly a photons would then escape through 
regions or patches with little neutral gas, and be blocked by gas 
at slightly higher column densities, while we measure our veloc- 
ities and covering fractions in the very densest clumps. 

For Haro 11, a picket-fence model of the ISM agrees 
well with the detection of the Lyman continuum (Ly C, < 
912A) escape dBerevall et al.l2006l:lLeitet et al.l201 ll). for which 
only direct paths are viable (see e.g. iHeckman et al.l [2001; 
IZastrow et al.ll201 IT) . An indirect method of estimating Ly C es- 
cape through the residual in the C n (A 1036) line points to that 
ESO 338 is also a strong candidate for Ly C escape (Leitet et al. 
2012, submitted), but unfortunately the redshift is too low for a 
direct measurement. 

The question still remains why, between knot B and C in 
Haro 1 1 with their relatively similar nebular dust content, Ly a 
is absorbed in knot B, which shows the larger kinematic offset. 
We hypothesize on the existence of diffuse remnant H i gas with 
a high covering fraction towards knot B. 
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Our measurements of the neutr al gas velocities in H aro 1 1 B 
agree rather well with the results in iKunth etaTI {l 998 ) where an 
outflow of ~60 km/s measured by low ionization state absorption 
lines in the UV was seen. However, the same spectra show Ly a 
in emission which would imply that some of the emission from 
knot C is included in the aperture as well. It is unfortunately 
not clear exactly where in th e galaxy the HST/ GHRS aperture 
was placed (see discussion in [Haves et alj |2007). Most likely, it 
was in between the knots, and it is quite like ly that the differen t 
components are mixed in the spectra. Indeed. lKunth et all (fl 998) 
see very broad absorption features (spanning roughly 200 km/s), 
indicating multiple absorption components. 

We note that a previous measurement of the neutral gas ve- 
locity in ESO 338 was made with HST/STIS in ISchwartz et al l 
(2006), which agrees well with our observations. The LIS ab- 
sorption features are weak and appear close to the systemic ve- 
locity. Unfortunately, the available HST/STIS spectrum for ESO 
338 is not optimized for exploring the Ly a emitting region that 
we investigate in this paper. In light of our results, it would be 
very interesting to explore both Haro 1 1 and ESO 338 with mul- 
tiple HST/COS and HST/STIS pointings to investigate the Ly 
a emission and absorption line profile in these galaxies in un- 
precedented spatial and spectral detail, which can be achieved 
by combining the two instruments. 

7. Conclusions 

We have presented the first spatially resolved measurements of 
the sodium doublet (Na D) in the two nearby Ly a emitting 
galaxies Haro 1 1 and ESO 338-IG04. Our results can be sum- 
marized as follows: 

- We find an outflow of neutral gas from knot B (which shows 
strong Ly a absorption) and slow infall towards knot C in 
Haro 1 1. In ESO 338-IG04 we find a slow or static interstel- 
lar medium. In the two latter cases, Ly a is seen in emission. 
The velocities that we find for the neutral gas are not what 
we would have expected from standard Ly a escape scenar- 
ios. Typically, a strong outflow is assumed to allow Ly a to 
escape more easily. 

- From the Na D line profiles we measure relatively high op- 
tical depths but small covering fractions of Na D (~ 10%). 
We estimate the minimum column density of Na D corre- 
sponding to our limits on the optical depth and find that Na 
D is likely only detected in the very densest clumps, where 
N(H i) > 10 20 cm -2 . Since Ly a is affected by resonant scat- 
tering already at column densities of N(H i) ~ 10 13 ~ 14 cm~ 2 , 
it is likely that the direct Ly a escape seen in e.g. knot C in 
Haro 1 1 and from ESO 338 is due to a picket-fence scenario 
where the interstellar medium is highly inhomogeneous and 
consists of both dense, neutral clumps as well as ionized gas 
along our lines of sight. 

- We see a larger kinematical offset in Haro 1 1 knot B than C, 
yet B shows strong Ly a absorption in contrast to the emis- 
sion from C. Given the relatively similar nebular dust content 
of the two knots, we hypothesize on the existence of a dif- 
fuse remnant H i component with a high covering fraction 
towards B, and a possible perpendicular outflow from C. 
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